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Two mesoscopic SQUID rings which are far from each other, are considered. A source of two- 
mode nonclassical microwaves irradiates the two rings with correlated photons. The Josephson 
currents are in this case quantum mechanical operators, and their expectation values with respect 
to the density matrix of the microwaves, yield the experimentally observed currents. Classically 
correlated (separable) and quantum mechanically correlated (entangled) microwaves are considered, 
and their effect on the Josephson currents is quantified. Results for two different examples that 
involve microwaves in number states and coherent states are derived. It is shown that the quantum 
statistics of the tunnelling electron pairs through the Josephson junctions in the two rings, are 
correlated. 
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I. INTRODUCTION 

Superconducting quantum interference devices 
(SQUIDs) exhibit quantum coherence at the macro- 
scopic leveli. This is a major research field within 
condensed matter, and has potential appHcations in the 
developing area of quantum information processing^''^. 
A lot of the work on superconducting rings investigates 
their interaction with classical electromagnetic fields. 

In the last twenty years nonclassical electromagnetic 
fields at low temperatures {k^T <C huj) have been studied 
extensively theoretically and experimentallj^ at both op- 
tical and microwave frequencies^ . They are carefully pre- 
pared in a particular quantum state, which is described 
mathematically with a density matrix p. The interaction 
of SQUID rings with nonclassical microwaves has been 
studied in the literatureSsi. In this case the full system, 
device and microwaves, is quantum mechanical and dis- 
plays interesting quantum behaviour. For example, the 
quantum noise in the nonclassical microwaves affects the 
Josephson currents. Experimental work, which involves 
the interaction of a Josephson device with a single pho- 
ton, has recently been reported*. 

An important feature of two-mode nonclassical mi- 
crowaves is entanglement. Entangled electromagnetic 
fields have been produced experimentally^. There is cur- 
rently a lot of work on the classification of correlated 
two-mode electromagnetic fields into classically corre- 
lated (separable) and quantum mechanically correlated 
(entangled)iS. In a previous publication— we have stud- 
ied the effects of entangled electromagnetic fields on dis- 
tant electron interference experiments. The interaction 
of entangled electromagnetic fields with two supercon- 
ducting charge qubits (that are approximated by two- 
level systems) has recently been studied^^. In that work 
it has been shown that the entanglement is transferred 
from the photons to the superconducting charge qubits. 
Related work has also been reported^'^. In this paper we 
study the effects of entangled electromagnetic fields on 
the Josephson currents of distant SQUID rings. 



We consider two mesoscopic SQUID rings, which are 
far from each other (Fig. 1). They are irradiated with 
entangled microwaves, produced by a single source. In 
this case the phase differences across the Josephson junc- 
tions are quantum mechanical operators. Consequently 
the quantum currents, which are sinusoidal functions of 
the phase differences, are also operators and their ex- 
pectation values with respect to the density matrix of 
the microwaves give the observed Josephson currents. It 
is shown that for entangled microwaves the currents in 
the two distant SQUID rings are correlated. We consider 
suitable examples of separable and entangled microwaves, 
which differ only by nondiagonal elements; and we show 
that the correlations between the induced Josephson cur- 
rents are sensitive to these nondiagonal elements. 

In Sec. II we consider a single SQUID ring and study 
its interaction with nonclassical microwaves. We assume 
the external field approximation, where the electromag- 
netic field created by the Josephson current (back re- 
action) is neglected. We also consider mesoscopic rings 
which are small in comparison to the wavelength of the 
microwaves. It is shown that under these assumptions 
the Josephson current is proportional to the imaginary 
part of the Weyl function of the nonclassical microwaves. 

In Sec. Ill we analyze the experiment depicted in Fig. 
1, where two distant SQUID rings are irradiated with 
entangled microwaves. We present examples of separa- 
ble and entangled microwaves that involve number states 
(Sec. IV) and coherent states (Sec. V). In Sec. VI we 
present numerical calculations for these examples. In Sec. 
VII we conclude with a discussion of our results. 



II. INTERACTION OF A SINGLE SQUID RING 
WITH NONCLASSICAL MICROWAVES 

In this section we consider a single SQUID ring and 
study its interaction with both classical and nonclassical 
microwaves. 



2 



A. Classical microwaves 

The current is /a — IiSihOa, where 9a — 2e^A is 
the phase difference across the junction due to the total 
flux $A through the ring. In the external field approx- 
imation, <i>A is simply the externally applied flux, while 
the back reaction (i.e., the flux induced by the SQUID 
ring current) is neglected. In other words the flux CI a, 
where C is the self-inductance of the ring, is assumed to 
be much smaller than the external flux <I>a. We consider 
a magnetic flux with a linear and a sinusoidal component: 

<^A^VAt + (j)A, (jjA^ Asin{ujit). (1) 

In this case the current is 

Ia = Ii sin[cjA^ + 2e^sin(wit)], ua = ScVa. (2) 

B. Nonclassical microwaves 

In this subsection we consider nonclassical electromag- 
netic fields, which are carefully prepared in a particular 
quantum state and are described by a density matrix p. 
In this case, not only the average values {E), (B) of the 
electric and magnetic fields are known, but also the stan- 
dard deviations AE, AB (and their higher moments) . 

A particular example is coherent versus squeezed mi- 
crowaves. In both cases the average values (E), (B) are 
sinusoidal functions of time, that can be made equal for 
suitable values of the parameters. However the uncer- 
tainties are AE — AB = 2^^/^ for coherent microwaves; 
and AE = a-^2-^^^,AB = a2-^^^, for squeezed mi- 
crowaves (where a is the squeezing parameter). 

Another way of describing nonclassical electromag- 
netic fields is through the photon counting distribution 
Pn — {N\p\N). For example, in the case of coherent and 
squeezed microwaves the distribution Pn is Poissonian 
and sub-Poissonian, correspondingly. One of our aims 
in this paper is to study how the quantum noise in the 
nonclassical fields, quantified by AE, AB, or with the 
distribution Pn, affects the Josephson currents. 

In quantized electromagnetic fields the vector potential 
Ai and the electric field Ei are dual quantum variables 
(operators). Strictly speaking the dual quantum vari- 
ables should be local quantities, but we consider meso- 
scopic SQUID rings which are much smaller than the 
wavelength of the microwaves. Therefore we can inte- 
grate these quantities over the SQUID ring and obtain 
the magnetic flux and the electromotive force: 

(j>= f Aidxi, Vemf = f Eidxi. (3) 
Jc Jc 

As explained above we work in the external fleld ap- 
proximation and we neglect the back reaction flux from 
the electron pairs on the external microwaves. In this 
case the flux operator evolves as 

(j){t) — —=[a) e-xp{iLot) + dexp(— iwt)], (4) 



where ^ is a parameter proportional to the area of the 
SQUID ring and the a\d are the photon creation and 
annihilation operators (e.g., Refsi^). 

In order to go beyond the external fleld approximation 
we need to consider the Hamiltonian 

i? = (^a^a + ^^ +i?SQUiD + ffi„t, (5) 

where -ffsQUiD is the SQUID Hamiltonian and Hu^t is 
the interaction term between the SQUID and the mi- 
crowaves. In this case the flux operator (t>'{t) evolves as 

= 4>{t) + --- (6) 

In this paper we work in the external fleld approximation 
and consider the flux operator of Eq. Q). 

Consequently the phase difference 9a is the operator 

^A = ^At + q[d^ e:sq){iujt) + aexp(— io^t)], q — \/2e^; (7) 

and the current also becomes an operator, 

Ia = Ii sinjciJA^ + q[d^ exp{iLut) + d exp(— ia;<)]}. (8) 

Expectation values of the current are calculated by tak- 
ing its trace with respect to the density matrix p, which 
describes the nonclassical electromagnetic flelds, 

(/a) = Tr(p/A) =/iIm[exp(zc^At)VF(AA)], (9) 
-^A = i(? exp(ici;ii). (10) 

Here W{x) is the Weyl functioni^ which is deflned in 
terms of the displacement operator D{x) as 

W{x) = Tt[pD{x)]- D{x) = exp(a;a^ - x*d). (11) 

The tilde in the notation of the Weyl function indicates 
that it is the two-dimensional Fourier transform of the 
Wigner function. 

In a similar way we can calculate the {l\) = Tr(/9/|). 
The second (and higher) moments of the current describe 
the quantum statistics of the electron pairs tunnelling 
through the Josephson junctions. As explained earlier, 
nonclassical electromagnetic fields are characterized by 
the photon counting distribution Pn = {N\p\N). The 
statistics of the photons threading the ring affects the 
statistics of the tunnelling electron pairs, which is quan- 
tified with the {I'a), {I a), etc. 



III. INTERACTION OF TWO DISTANT SQUID 
RINGS WITH ENTANGLED MICROWAVES 

In this section we consider two SQUID rings far apart 
from each other, which we refer to as A and B (Fig. 1). 
They are irradiated with microwaves which are produced 
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We also calculate the second moments 



/2Tr[pA(sin^A)'], 
(/|) = /|Tr[pB(sin0B)']. 



(18) 
(19) 

As explained earlier, the statistics of the photons thread- 
ing the ring affects the statistics of the tunnelling electron 
pairs, which is quantified with the (/a/b), (-^1)' (-^b)' ^t*^- 
In the following sections we consider particular ex- 
amples for the density matrix p of the nonclassical mi- 
crowaves that interact with the two SQUID rings, and 
examine its effect on these quantities. 



IV. MICROWAVES IN NUMBER STATES 



FIG. 1: Two distant mesoscopic SQUID rings A and B are 
irradiated with nonclassical microwaves of frequencies uoi and 
UJ2, correspondingly. The microwaves are produced by the 
source Sem and are correlated. Classical magnetic fluxes Vxt 
and Vei are also threading the two rings A and B, correspond- 
ingly. 



by the same source and are correlated. Let p be the 
density matrix of the microwaves, and 



PA = TrBP, Pb = TrAP, 



(12) 



the density matrices of the microwaves interacting with 
the two SQUID rings A, B, correspondingly. When the 
density matrix p is factorizable as pfact — PA ® Pb the 
two modes are not correlated. If it can be written as 
Pscp = ^iPiPM ® PBi, where Pi are probabilities, it is 
called separable and the two modes are classically corre- 
lated. Density matrices which cannot be written in one 
of these two forms are entangled (quantum mechanically 
correlated). There has been a lot of work on criteria 
which distinguish separable and entangled statesi^. 
The currents in the two SQUIDs are 



(/a) = /iTr(pAsin6'A), 
(/b) - /2Tr(pBsin^B). 



(13) 
(14) 



The (/a) is written in terms of the Weyl function W{Xa) 
in Eq. lO, and similarly for B one may obtain (/b) = 
/2lni[exp(zcjBi)II^(AB)], where Ab = iq e:sq){iu!2t) and 
ujb = ScVb. 

The expectation value of the product of the two current 
operators is given by: 



(/a/b) = /i/2Tr(psin0Asin0B). 
We consider the ratio of the currents 

(/a/b) 



R = 



(/a)(/b) 



(15) 



(16) 



For factorizable density matrices pfact = Pa®Pb we easily 
see that /?fact — 1- For separable density matrices psop = 
J2iPiPAi PBi we get 



^scp 



J2iP'i{^Ai){lBi 



(EfePfc(/Afe))(E/K(/B()) 



(17) 



Pscp = T;iWiN2){NiN2\ + \N2Ni){N2Ni\), (20) 



We consider microwaves in the separable (mixed) state 
1 

where A^i ^ N2. We also consider microwaves in the 
entangled state \s) = 2-^/'^{\NiN2) + \N2N1)), which is 
a pure state. The density matrix of |s) is 

Pent = Pscp + \{\NiN2){N2Ni\ + \N2Ni){NiN2\), (21) 

where the Psop is given by Eq. (|20|l . It is seen that 
the Pent and the Psep differ only by the above nondiago- 
nal elements, and below we calculate their effect on the 
Josephson currents. 

We note that it is possible to have 'interpolating' den- 
sity matrices of the form 



Pp 



PPscp + (1 - J3)Pont 



(22) 



^ Pscp 



1-p 



{\NiN2){N2Ni\ + \N2Ni){NiN2\) 



where < p < 1. Below we present results for the two 
extreme cases of Psop, where the nondiagonal terms make 
no contribution; and for the Pent, where the nondiagonal 
terms make maximal contribution. We also present nu- 
merical results for the case of pp. 

In this example, the reduced density matrices are the 
same for both the separable and entangled states: 



Psop, A — Pent, A ~ Psep,B — Pent,B 

= l:m){N,\ + \N2){N2\). 



(23) 



Consequently in this example (/A)sep ~ (/A)ent, and also 



'B/sep 



'B/ent- 



A. Classically correlated photons 

For the density matrix psep of Eq. H2U|I we find 

(/a) - /iC'sin(wAt); (24) 
(/b) = hCsiniujBt); (25) 

C = lc^v(-^)[LNAq^) + LNAq% (26) 
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where the i"(a;) are Laguerre polynomials (in the case 
of Eq. we have a = 0). The currents (/a), (/b) 

are in this case independent of the microwave frequencies 

The expectation value of the product of the two cur- 
rents [Eq. lfT3|l ] is 

(/a/b)sop = /i/2Cisin(wAt)sin(cjBi), (27) 
Ci = expi~q^)LNAq^)LN,{q^). (28) 

Consequently the ratio R of Eq. (|16|) is 



^ Ci ^ 4£Ar,(g^)£jv,(g^) 

''^''P C2 [L^,(q2) + LAr,(92)]2- 



(29) 



In this example the i?scp is time-independent. 

The moments of the currents, defined by Eqs. (|18|l . 
((T^ . are also calculated: 

(/i) = :|[l-C2Cos(2c^Ai)], (30) 

(/i) = :|[l-C2Cos(2c^Bi)], (31) 

C2 = ^eM~2q^)[LN^iM^)+LNAM% (32) 



B. Quantum mechanically correlated photons 

For the case of pcnt the reduced density matrices pA, Pb 
are those given by Eq. H23|l . and consequently the 
(/a),(/b) are the same as in Eqs. (1^ . and the 

(il), (/|) are the same as in Eqs. ^U^, (1^ . 

However in this case the (/a^b) is 



(4/1 



B/cnt 



(/a/b>s 



(33) 



where 



•OSS = -hhCs COs[(iVi - N2)iuJi - UJ2)t] 
X [cOs(ti>Ai + ^Bt) — ( — l)^i~^2 cos(ti>Ai — ^Bt)], 



N2 



(34) 
(35) 



The term /cross is induced by the nondiagonal elements 
of Pent J and depends on the photon frequencies U!i,uj2- 
This term quantifies the difference between the effect of 
separable and entangled microwaves on the Josephson 
currents. We note that the nondiagonal terms of pont [Eq. 
(|21|l ] are small and consequently they are very sensitive 
to the back reaction. Therefore our results which neglect 
the back reaction are relevant to experiments with small 
Josephson currents. In other words it is required that 
the fluxes £a/a and >Cb/b are much smaller than the 
external flux. 

The ratio R of Eq. Hlt)|) can be simplified in two 
distinct expressions according to whether the difference 



Ni — N2 is even or odd. In the case Ni — N2 = 2k, the 
ratio is 

(2fe) 4L^2fe(^2)^2fc (^2) 

/?e„t =^-P+[i^^(;2)+i^^\^2)]2C°«(^^)' (36) 

where 

n = {Ni-N2){oJi-u}2)- (37) 

It is seen that the /?c„t oscillates around the i?scp with 
frequency fl given by Eq. 1)3 7() . If there is no detuning 
between the nonclassical electromagnetic fields, i.e. lui — 
UJ2, then /?c^t' is constant, although it is still Rem 7^ Rscp- 
In the case TVi — iV2 = 2fc -f 1 the ratio is 

p(2fc+l) _ p 
-"-cnt ~ -n,sop 



iL-^f-Hq')Ll\+H<l') cos(m) 
[Ljvi (g^) + Ln^ (q2)]2 tan(tJAi) tan(cjBt) 



• (38) 



In both cases the /?ont is time-dependent and it is a func- 
tion of the photon frequencies a;i,a;2, in contrast to the 
case of /?scp (which is time- independent) . 



V. MICROWAVES IN COHERENT STATES 

We consider microwaves in the classically correlated 

state 



p,ep - -(IA1A2XA1A2I + 1^12^1) (^2^1 1). (39) 



The 1^1), 1^2) are microwave coherent states (eigenstates 
of the annihilation operators). We also consider the en- 
tangled state \u) = J\f{\AiA2) + [^2^1)), with density 
matrix 

Pont = 2A^^Pscp 

+ M^i\A,A2){A2A,\ + \A2A,){A,A2\), (40) 
where the normalization constant is given by 

AA= [2 + 2exp(-|.4i- A2|2)]"'/'. 



(41) 



A. Classically correlated photons 

For microwaves in the separable state of Eq. (|39|l the 
reduced density matrices are 



Psep,A = Psop,B = 2(1^1X^11 + \^2)(A2\), (42) 

and hence the currents in A and B are 

J 2 

(/a)scp = Y'^'^P( — ^){sm[ujAt + 2q\Ai\cos{ujit- 9i)] 



sin[a;At + 2(7|^2| cos(a.)ii - 6*2)]}, 



(43) 



5 



(-fB)scp = yf^xpl — ^){sm[ujBt + 2q\Ai\ cos{uJ2t - 9i)] 
+ sm[LUBt + 2q\A2\cos{LU2t-e2)]}, (44) 

where 6i = arg(Ai), and 62 = arg(A2). We have also 
calculated numerically the ratio i?scp- 

B. Quantum mechanically correlated photons 

For microwaves in the entangled state of Eq. H4()|l the 
reduced density matrices are 

/Ocnt.A = Pent,B^J^^{\Ai){Ai\ + \A2){A2\ 

+ t\Ai){A2\+t*\A2){A,\), (45) 

where 

r = {AM2)=cxp(^-\^-\^+AlA2y (46) 
The current in A is 
(/A)c„t = 2M^{iA)scp+M^EFi exp (^-^^ h, (47) 
where 

E = exp[-|^i|2 - + 21^1^2! cos(e'i - 6*2)], (48) 

Fi = {exp[g|Ai|5^,i(i) ~ q\A2\SAMt)] 
+ eM~q\Ai\SAAt) + QlMSAAt)]} 
X siniujAt + q\Ai\CAAt) + qlMCAAt)], (49) 

and 

Saa = sin(wii - 6*1), S'a,2 = sin([Jii - 6*2), 
Ca,i ~ cos{iLJit — 9i), Ca,2 = cos{iLjit — 92)- (50) 

The current in B is 

(/B)cnt = 2N'^{iB)scp+M^EF2exp (^-^^ I2, (51) 

where 

F2 = {exp[g|Ai|5B,i(i) - <z|^2|5s,2(t)] 
+ exp[-g|Ai|^B,i(0+(7|A2|5s,2(t)]} 
X sin[wB< + g|^i|Cs,i(t) + q\A2\CBAt)h (52) 

and 

Sba = sin(w2t - 6*1), Sb,2 ^ sm{LU2t - 92), 
Cba = cos{ijj2t - 61), C_B, 2 = cos(w2t - 6*2). (53) 

We have also calculated numerically the ratio i?ont- 

As we already explained, the nondiagonal terms in pont 
are very sensitive to back reaction and therefore these 
results are relevant to experiments with small Josephson 
currents. In other words it is required that the fluxes 
CaIa and CbIb sue much smaller than the external flux. 



VI. NUMERICAL RESULTS 

In all numerical results of Figs. 2 to 6 the microwave 
frequencies are wi = 1.2x10"'*, 0^2 = 10"*, in units where 
kB = h = c = 1. The critical currents are Ii = I2 = 1- 
The other parameters are ^ = 1, lua — '^-'i, 'j-'b — '-^2, 
Ni = 2, N2 = 0, and the arguments of the coherent 
eigenstates are 9i — 92 = 0. For a meaningful comparison 
between microwaves in number states and microwaves in 
coherent states, we take = iVi and |^2p = N2, 

so that the average number of photons in the coherent 
states is equal to the number of photons in the number 
states. 

In Fig. 2 we plot i?sep against {uji — 0^2 )t for currents 
induced by microwaves in the number state of Eq. H20(l 
with A^i = 2, Af2 = (line of circles), and the coherent 
state of Eq. with Ai = ^2,^2 = (solid line). It is 
seen that two different microwave states with the same 
average number of photons give different results on the 
quantum statistics of the electron pairs. 

In Fig. 3 we plot i?sep — -Rent against (uji — 1^2)1 
for currents induced by microwaves in (a) the number 
states of Eqs. ^ and (EB with A^i = 2, iV2 = 0, 
and (b) the coherent states of Eqs. and (|in|l with 

Ai — ^/2,A2 = 0. It is seen that the separable and 
entangled states, which differ only by nondiagonal ele- 
ments, give different results. As expected, the difference 
(which shows the effect of the nondiagonal elements) is 
small, but it is nonzero. 

In Fig. 4 we plot (a) (/a)sop - (-fA)ent, and (b) 
(-^a)scp — (-^A)ent: against {uji ~uj2)t for microwaves in the 
coherent state pA.scp of Eq. and pA,cnt of Eq. 
with Ai = a/2, ^2 — 0. For coherent states Pcnt.A is not 
equal to pscp.A [cf. Eqs. (|^ . (021)] and consequently the 
corresponding currents are different. For number states, 
the currents corresponding to Pcnt,A and Pscp,A are the 
same because Pcnt.A — Pscp.A [cf. Eq. H23|l ] . 

In Fig. 5 we plot (/A-fB)sop - (-^A-?^B)cnt for (a) number 
states of Eqs. and ^ with iVi = 2, iV2 = 0, and 

(b) coherent states of Eqs. lIH^ and with Ai = 
\/2, A2 — 0, against {ijji — uj2)t. In this figure also, we get 
different results due to the nondiagonal elements in the 
entangled state. 

In Fig. 6 we plot the ratio Rp of Eq. l(T?)|l against p e 
[0, 1] for currents that are induced by the interpolating 
density matrix p-p of Eq. H22I) for number states with 
A^i = 2, iV2 = 0. The time t has been fixed so that 
(wi - U2)t = 7r/2. 



VII. DISCUSSION 

We have considered the interaction of SQUID rings 
with nonclassical microwaves. We have assumed the ex- 
ternal field approximation, where the electromagnetic 
field created by the Josephson currents (back reaction) 
is neglected. We have also considered small rings in com- 
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parison to the wavelength of the microwaves and taken 
as dual quantum variables the magnetic flux and electro- 
motive force of Eq. Q . 

The Josephson current is an operator and its expecta- 
tion value with respect to the density matrix of the non- 
classical microwaves is the observed current. We have 
shown that the expectation value of the current is pro- 
portional to the imaginary part of the Weyl function [Eq. 
©]. This shows clearly how the full density matrix of 
the microwaves affects the Josephson current. The higher 
moments of the current {I^) can also be calculated and 
used to quantify the statistics of the tunnelling electron 
pairs. It has been shown that the statistics of the irradi- 
ating photons determine the tunnelling statistics of the 
electron pairs. 

We have also considered the interaction of two distant 
SQUID rings A and B with two-mode nonclassical mi- 
crowaves, which are produced by the same source. It has 
been shown that classically correlated (separable) and 
quantum mechanically correlated (entangled) photons in- 
duce different Josephson currents and different tunnelling 



statistics in the two devices. The results show that the 
entangled photons produce entangled Josephson currents 
in the distant SQUID rings. This can have applications 
in the general area of quantum information processing. 

The work can be extended in various directions. The 
first is to take into account the back reaction and include 
the extra terms which we have neglected in Eq. © . This 
can be done numerically. The second direction is the 
study of Bell-like inequalities for the Josephson currents; 
which are violated when the currents are entangled. An- 
other direction is the potential use of the system as a 
detector of entangled photons. There is a lot of work on 
the use of mesoscopic devices as detectorai^. Applica- 
tion of the present work in this direction requires further 
work, which could lead to the development of a detec- 
tion system for entangled photons based on two distant 
SQUID rings. 



* Corresponding author: d.i.tsomokos@physics.org 
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FIG. 2: -Rscp against {loi — uJ2)t for the number state of Eq. 
I|2()|l with A''i = 2,N2 = (Une of circles), and the coherent 
state of Eq. ^ with = ^2,^2 = (solid line). The 
photon frequencies are uji — 1.2 x 10~* and UJ2 — 10~*, in 
units where fcs = h = c = 1. 




FIG. 3: -Rscp — Rcta against {u)\ — uj2)t for (a) number states 
of Eqs. ((213 and ^ with iVi = 2, N2 ^ 0, and (b) coherent 
states of Eqs. ^ and ^ with Ai = V2,A2 ^ 0. The 
photon frequencies are ui — 1.2 x 10"* and L02 ~ 10"'*, in 
units where ks ~ h — c — 1. 
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FIG. 4: (a) (/a)scp - (/a)™!, and (b) {/l)scp - {/l)ont against 
{uji —uo2)t for the coherent state psop.A of Eq. H42|l and pcnt.A 
of Eq. 1451 with A\ = \/2, A2 = 0. The photon frequencies 
are uji = 1.2 x 10"'* and uj2 = lO"**, in units where ks ~ h = 
c = 1. 
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FIG. 5: (/a/b)scp — {/A-?B)cnt for (a) number states of Eqs. 
iCT and ||2IJ with TVi = 2, iV2 = 0, and (b) coherent states of 
Eqs. and igHJ with Ai = V2,A2 = 0, against {uji —U2)t. 
The photon frequencies are uji = 1.2 x 10~* and ijJ2 ~ 10~*, 
in units where kB = h — c = 1. 
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FIG. 6: Rp against p for currents that are induced by the 
interpolating density matrix pp of Eq. H22|l for number states 
with A'^i —2, N2 = 0. The time t has been fixed so that 
[uji —uj2)t — -K 12. The photon frequencies are uj\ = 1.2 x 10~* 
and (jj2 = 10^*, in units where fcs = ?!, = c = 1. 



